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ABSTRACT: The following report io concerned with a method for 

the approxlxnato calculation of compressible flows 
about profllos with local regions of supersonic 
velocity. The flow around a slender profile Is 
treated as an example . 
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velocity. 
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I. STATEMENT OF THE PROBLEM 


Several methods are known, for the calculation of compressible 
flows at high subsonic velocities. The resulting approximate 
solutions are quite useful as long as sound velocity Is not 
exceeded at any point of the flow field. However, apparently all 
these approximate calculations, i/ithout exception, cease to 
converge or to render useful flow patterns If the condition of 
purely subsonic flow is no longer satisfied. Moreover, numerous 
tests In wind tunnels confirmed the result that the reconversion 

*”Berechnung kcmpresslbler Strctaungen mlt b'rtlichen TTber- 
schallf eldem. " Zentrale f^r wissenechaftliches Berlchtswesen 
der Luftfahrtforschung des Generalluftzeugmei stars (ZWB) Berlln- 
Adlershof , Forschungsberlcht Nr. 1794, Berlin- Adlershof, dan 
7. August 1943. 


e 


TtACA. TM Ho. 1114 


of local super sanlo flaws Into subsonic flows la real flows 
takas place not steadily but generally by means of a compression 
shock which completely changes the whole flow pattern. 1 There- 
fore, the question arises whether a continuous process from 
supersonic to subsonic flow Is at all possible for bodies In 
parallel flow, even though the formation of a boundary layer on 
the surface of the body la at first neglected. 

An approximation method was tested In the OTL In order to 
clarify these questions; this method makes possible the calculation 
of flows with local regions of supersonic velocity. 


II. SURVEY OF 2BB METHOD USED 


The development of the method started from the fact that the 
known approximation methods for compressible flows, which are 
without exo option based cm a step-by-stop improvement of the 
Incompressible flow, are quite appropriate for the limited domain 
of the pure subsonic flow extending from Infinity to the sonic 
velocity boundary In the flow field near the body, however, 
for the region of the local supersonic flow a method based on 
the properties of supersonic flow will be subsequently used, as 
for instance the method of characteristic s of Prandtl-Buaemam. 
Accordingly the partial areas of supersonic and subsonic flow, 
respectively, are calculated separately by different methods 
which in each case are treated according to the peculiarities of 
the partial flow to be calculated. The partial flows that were 
thus determined must them 'bo Joined In such a manner that the 
flows agree on the surface of contact of the two regions, that Is, 
on the sonic velocity boundary, with respect to magnitude as well 
as to direction of the velocity. Figure 1 shows a schematic 
representation of the boundary between the two flow areas. 


1. Subsonic Region with 3na.ll Perturbation Velocities 

Prandtl 1 s rule will represent a good approximation for a 
great part of the outer subsonic flow area, since the perturbation 
velocities caused by the profile are sufficiently small up to 
seme distance from the sonic boundary so that the assumptions of 

account will be taken here of flows which slightly exceed 
sound velocity because for them there is no certain, way cf 
distinguishing between a steady transition and a compression 
shock; neither by pressure -distribution measurements nor by optical 
observations, for Instance according to the Bchlleren method. . 
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Prandtl’ a approximation are rather well satisfied. The subsonic 
flov of this region can, of course, not he obtained by simply- 
distorting the incompressible flov around the profile according 
to Prandtl* s rule. First, the local supersonic velocity field 
requires essentially more space than the subsonic flov. The 
outer subsonic streamlines, therefore, are widened outward not 
only by the body In the flov but also by the additional displace- 
ment due to the supersonic velocity field. However, this additional 
widening of the streamlines by the supersonic velocity field can 
only affect the subsonic flov like a modification of the boundary 
conditions at the Bonic velocity boundary; such a modification 
may be represented In a simple way by source end sink distributions, 
dipoles, end so forth. 


2. Subsonic Region with High Porturbcvtion Velocities 

As mentioned before, the part cf the mb sonic flow which has 
to be calculated according to Prandtl does net extend as far as 
the sonic boundary; it only reaches up to a boundary lino near 
the sonic boundary determined by an agreed sufficiently small 
value of the perturbation velocities. (Tills boundary line is 
represented in fig. 1 by a dashed lino.) Tho calculation for 
the region from this line to the sonic velocity boundary must 
generally be carried out by an improved subsonic method. A 
numerical method of calculation seems to he partici’larly appro- 
priate for this intermediate region. Tho outer flow according 
to Prandtl, which io assumed as known, may by this method by 
continued a little further, namely up to the sonic velocity 
boundary. 

One assumes, for instance, that the flow field is covered by 
a rectangular grid of selected points. The velocity component s 
v x and Vy at these points in the outer field are known for 
each case according to. Prandtl* s approximation . Then the 
velocities at the inner grid points may be calculated from the 
known values at the outer grid points if the differential 
quotients which are decisive for this continuation are appro- 
ximated by the corresponding difference quotients. 2 Figure 2 
shall be considered as an example. Magnitude and direction of 
the velocities at all grid points outside of the boundary lines 
are assumed as known. The air density shall be known also. 

Then the exact equations of the continuity and the irrotatianality 
are, for the indicated grid point "2": 

2 I wish to express my thanks here to Dr. H. Schubert /DVL 
for hia suggestion to treat compressible flows according to 
the difference method. 
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or, In the notation of difference quotients: 

r x x — — ^L. + p Tj A Vi 4 + v T x P 3 **1 + p 0 x 7 3 y l a o 
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If one assumes the grid to he quadratic with meshes of equal width 
then As n Ay and the equations given above are simplified as 
follows: 


t i 2 x %,U + P 2 x *Va,4 - -Vs * *\l ' p 2 x A V 3 ,1 (la) 

*Va,4 ’ At *3,1 (») 

From the second equation one Immediately obtained the velocity 
component v yA at the point A, while the first equation gives 

a relation between the velocity component v and the air 

density at the point A. Velocity and air density are connected 
in a rather complicated way by the adiabatic relation 



(p Q hefe represents the air density in the gas at rest a* the 
critical sound velocity) . Therefore, it will be practical to 
use the linearized relation Ap ■ ^ Av instead of the exact 
function., the error being negligible, beoause of the small width 
of the meshes: 
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Compare figure 3. 
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The Telocity ooinpanents v and v_ for the new grid point A 

x A . „A 

ero known, according to those calculations, from equations (1 b) and ( 3) . 
Therefor©- the air density also is known according to equation (2). 

Thus the intermediate region "between tho subsonic flow with mjn.1 .1 
perturbation Telocitios (calculation according to Pranitl) and 
the houndary may he determined and the position of the sonic 
Telocity boundary and the Telocity direction along that boundary 
will he obtained as final result. It will have to ho further 
invest! gated whether a calculation of the whole supersonic velocity 
field beyond the sonic boundary up to the surface of the "body by 
this method would be practical. In the main, the time spent on 
the calculations and the accuracy of the method which has been 
developed here as contrasted with the method of characteristics 
of Prandtl -B ub stnaim will have to be considered. 

The difference method which has been developed here does not, 
howerer, permit a start of the numerical calculation at an arbitrary- 
distance frcm the body in the flow; for at a very great distance 
from the body the calculation would not be sufficiently accurate 
since tho meshes of the grid would be too small. This conclusion 
can bo drawn from the fact that at a great distance the flow about 
any body without external forces can be replaced, except for email 
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deviations, "by the flow field around a single dipole. It Is 
precisely the mail deviations from the dipole flow, which exist 
even at the greatest distance, that permit a continuation of the 
flow, for Instance according to the difference method, not about 
a cylinder hut around the special profile. 


3. gupereonic Region 

From the calculations described above, which were carried 
out In the subsonic area according to Prandtl , s rule, or by the 
difference method, the boundary line of the supersonic velocity 
field was obtained. The velocity along this lino equals the 
sound velocity and the direction ef the flow at every point of 
this line Is known. The next problem consists In Joining the 
corresponding supersonic velocity field to this boundary line; 
for instance, by the method of characteristics of Prandtl- 
Busemann. According to the method of characteristic a It is 
known that the changes of state occurring In a supersonic flow 
are manifested In expansion cr compression waves. Therefore, 
when the flow passes through these lines of disturbance, the 
flow velocity and the flow direction are modified by definite 
values which can easily be determined from the graphical 
representation of the characteristics. 3 It can be Inferred from 
the condition of constant velocity at the sonic boundary line 
that two waves must originate at every point of the sonic ^pundary, 
namely an expansion and a compression wave. (See fig, b.) The 
velocity direction Jumps from one starting point of these waves 
at the sonic boundary to the next by a fixed amount for each case. 
Since these starting points of the waves can be shifted in any 
way along the sonic velocity boundary, the beginning of a 
graphical representation of characteristics may be drawn for 
any given distribution of direction on the sonic boundary. 

Likewise the starting points of the waves may be shifted not 
only on the sonic boundary but also towards points outside of 
the sonic boundary, so that the graphical representation of 
characteristics may be adjusted to any shape of the sonic 
boundary line. Therewith It has been demonstrated that at least 
in certain cases there exists a possibility of continuing the 
subsonic flow into a corresponding supersonic flow; finally there 
results the new profile contour of the body in the flow,' 1 

3 Ccmpare L. Prandtl: "Fuhrer durch die Stramungslehre . " 
Verlag Vioveg and Sohn, Braunschweig (1942) S. 257 und folgende. 

*In order to render possible sufficiently accurate drawings 
of grids of characteristics for local supersonic velocity field®, 
a graphical representation of characteristics with an interval 
of l/5° has been completed at the BVL; copies may be obtained. 
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IH. CALCT3LATI0B OP AB EXAMPLE POE A MIXED 
SUPERSONIC -SUBSONIC PLOW 


A simple example will demonstrate the combined effect of the 
calculation methods for various flow regions. The velocity field 
for a two-dimensional slender tody was determined In the proximity 
of the "body by means cf conformal transformation and the associated 
distortion according to PrraoAtl for the Mach number M a 0.86.- 
Thus there resulted a sonic velocity boundary as indicated in 
figure 5(a)- This eanlo boundary, which was obtained according 
to Prandtl* s method, could be further Improved by continuing the 
flow step by step to th% Improved sonic velocity boundary; one 
would have to start from a boundary lino with sufficiently small, 
perturbation velocities using the difference method described 
above. This correction of the sonic velocity boundary was at 
first disregarded for the sake of simplicity, since in general, 
it can bo neglected for the basic calculations plenned for this 
report. By means of the directions of velocity along the sonic 
boundary, the neb of characteristics for the supersonic flow can 
be drawn so that the contour of the body In the flow Is a stream- 
line of this supersonic field. Only expansion waves start from 
the surface of the body, while compression vaves, which are 
necessary for a reconversion of the supersonic flew Into a sub- 
sonic flow start without exception from the sonic boundary. This 
behaviour Is an e scent ial criterion of the grid of characteristics. 
The manner of reconversion of the supersonic into a subsonic flow 
(that is, for Instance, whether steadily or by means of a cam- 
pro a a Ion shock) Is decisively influenced by the sonic velocity 
boundary. 

Of course, the contour of the body at the boundary of the 
supersonic field which has been determined In the foregoing 
calculation of the example can not coincide with the contour of 
the original body which was originally distorted according to 
Prandtl, since the supersonic field, because of the expansion 
of the air, requires more space than the flow field which was 
distorted according to Prandtl. Therefore the surface of the 
body in the Buper sonic region will be flatter than the surface 
of the original body In order to satisfy this increased need for 
space. The outer subsonic flow must be adjusted to this Increased 
need of space of the supersonic field In order to obtain the 
original oontour as the result of the calculation. This adjustment 
may be made by replacing the supersonic field with dipoles or 
source-sink distributions, which will widen the streamlines of 
the subsonic region originally obtained by Prandtl* s method. The 
strength of this source-sink distribution Is known from the 
condition that the space which Is to be added has to cover tho 
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difference of the flov density In euperaonio flov and flov 
according to Prandtl and that this space corresponds In 
magnitude to the flattening of the original "body which resulted 
In figure 5(a). By a stop-hy-step approximation the contour 
which results at the end of the calculation may thus he 
adjusted to the desired ocntour. 

However, the example treated here already permits the 
conclusion that vhon the local Bound velocity is exceeded there 
exist solutions with reconversion of the local super sonic 
velocities Into subsonic velocities without compression shock, 
at least In frlctlonlees flow. 

The pressure distribution which was fo’md for the example 
treated. - Is compared to the pressure distribution of another 
slender body (fig. 5(bJ vhloh was determined according to the 
same law of conformal mapping and should agreo woll with the 
body of the example. (Sue fig. 5(&)-) The extrv^ne flatness 
of the pressure distribution In the supersonic region is a 
special characteristic for the body with local region of sivper- 
Eonic velocity. This behaviour corresponds to the strong 
curvature of the contour of the body Immediately after the 
flow has entered the supersonic region. Jhio to this large 
local curvature, the body B In this region will produce even 
■for incompressible flow higher negative prop sure a than the 
body A. For this reason the two pressure distributions which 
were graphically represented cannot be directly compared. 


TV. APERCJX3MAEE TREA3MKRT OF LOCAL EEGIONS 
OF SUPERSONIC VELOCITY 


The representations of characteristic a considered thus far 
were notable for the fact that expansion or compression waves 
running In the same direction never Intersect. (Compsro fig. 5(a).) 
An Intersecting of waves running In the same direction Is basically 
impossible for expansion waves since such waves always diverge 
from their starting point. But compression waves convorga; there- 
fore^ they could well form an envelope and cause a compression 
shock. 

Such overlappings of waves were found for the sample body of 
figure 5(a) when the free- stream velocity was increased fr o m 
M a 0.86 to M «» 0.90. The conotruction of the expansion waves 
was started at the sonic boundary (In opposition to their actual 
direction) and already for the Mach number M b 0,85 (according 
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to figure 5(a)) led to a pressing together of the expansion, 
waves near ' the "contour; fotf the Mach number M'b 0^90 they 
ran Into each other even "before reaching the contour. This 
condition Is physically senseless, that Is, a symmetrical 
solution for the sonic "boundary which was found 'according to 
the Indi cated scheme is, at this" Mach number M ■ 0„90, no 
longer possible. In this case the "expansion chock, " according 
to a consideration of Prandtl, rather will dissolve into a 
group of divergent expansion waves originating at the contour, 
while the c ompre ssion shock which limits the supersonic area 
persists. In this maimer one easily obtains the flow picture 
found from many tests for. which the local supersonic flow is 
no longer reconverted into the subsonic flow through a continuous 
phenomenon but by means of a compression shook. 

The preceding treatment was based upon the sonic velocity 
boundary line which was obtained according to Prand5J. , s rule 
without additional source- sink todies which woi’la have corre- 
sponded to the greater need for space of the supersonic region. 
The resulting solution which was physically senseless might 
well be caused by the perhaps no longer appropriate conditions 
for the existence of the sonic velocity boundary line. In order 
to obtain a clearer view one considers the limiting case where 
the expansion waves converge exactly at the contour . A flow 
around a comer develops, as represented schematically in 
figure 6; In the supersonic region such a flow is pooslble 
without flow separations. Thin assumption and a given sonic 
velocity boundary offer the supersonic region the opportunity 
to fill the largest possible space without concave curvature of 
the body in the supersonic domain. 

An admission of ooncave surfaces would presuppose a sonic 
velocity boundary which could no longer be produced from an 
Incompressible flow about a profile without additional source 
and sink distributions In place of the supersonic region; or 
else one would have to drop the condition of flow symmetry. 

Figure 5(a) demonstrates clearly how the surface always shows 
a convex curvature when the Incoming wave is a compression wave, 
the outgoing wave an expansion wave. However, for the opposite 
condition of concave curvature In symmetrical flow the coming 
In of an expansion wave and going out of a compression wave 
would result In the amission of expansion waves from the sonic 
boundary In both cases. But an expansion wave which starts at 
the sonic boundary cannot be obtained by simply exchanging the 
compression and expansion waves, for Instance according to the 
diagram In figure h; for there must always first appear an 
expansion wave and than a compression wave in tho direction of 
the flow at the sonic velocity boundary, beoause the flow can 
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pass from the sonic ■boundary into the supersonic region only 
by means of an expansion. An expansion wave starting at the 
sonic boundary is possible only for a certain shape of the 
comers in the graphical representation of characteristics: 
two expansion waves are sent out from the comer simultaneously. 

( Comp are fig. 7.) There is, however, a condition: the sonic 

boundary In this region must take a steeper course than the 
expansion wave which is inclined toward the direction of the 
flow by the Mach angle. The streamline through such a comer 
point lies between the expansion waves, while for a noimal comer 
it is tangent to the apex of the triangle formed by the expansion 
and compression wave. 

It la an Important criterion for the comer with two expansion 
waves that the streamlines become steeper, not flatter, with 
Increasing distance from the body; otherwise the flow would not 
fit together after transition throught the two expansion waves. 

This result agrees with the fact that In a symmetrical supersonic 
region also the direction of the streamlines grows steeper with 
Increasing distance from the profile. The cause of this phenomenon 
Is that the streamlines (because of the mATrirmnn flow density at 
sonic velooity) are closest to each other at the sonic boundary 
while there will be the greatest distance between them for the 
supersonic region at the point of maximum velocity, that is, 
generally, at the point of greatest thickness of the body. 

Such distributions of slopes occur only In supersonic regions; 
they are not possible In subsonic areas for a flow about bodies 
with convex contours. Therefore it is net surprising that the 
sonic velocity boundary, which was obtained by distortion of the 
Incompressible flow about the body, will not lead to a useful 
supersonic flow pattern as long as there are no singularities, 
as for Instance source-sink distributions in the Incompressible 
flow by means of which the effeot of the Increased, expansion of 
the air In the supersonic region can be calculated. 

Presently a method is under investigation in which the sonic 
velocity boundary may be adjusted to tho increased ncei for 
space of the supersonic flow by arrangement of singularities in 
the flow; the results of these Investigations will soon be 
published separately. 


V. SOWETRICAL AND ASYMMETRICAL REGIONS OP SUPERSONIC PLOW 


As stated In detail In the preceding chapter, concave body 
surfaces In the supersonic region are dependent upon the particular 
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typo of comers In the net of characteristics ■which Is characterised 
by two expansion waves starting from a point of the sonic boundary. 
(Camp ore fig. 7.) The stipulation of* flow syinaetry then causes 
compression wares to start from the corresponding points of the 
contour which in turn cause these compression wares to proceed 
toward the sonic boundary in pairs . Since these compression wares 
which start at the contour always converge, such solutions ere 
possible only as long as the radiated compression wares do not 
Intersect; otherwise a compression shock will develop. Therefore 
the existence of symmetrical flows with local supersonic regions 
cannot be oounted upon for strongly concave curvatures or for high 
Mach numbers, because the compression shock would presuppose a 
symmetrical expansion shock In a corresponding location; this 
expansion shock, however, is not physically possible. 

Except for the special case described above, supersonic flows 
an principle tend toward flow symmetry as will be shown in the 
example represented in figure 8(a). There will always develop a 
symmetrical flow along the wall with two changes in direction for 
incompressible flow without separations. It is, however, known 
that for pure supersonic flow the flow along the ■trail is unsynmetriceDL 
(fig. 8(a)), if no additional guiding surfaces exist. Conditions 
of symmetry can be achieved through, calculations for the mixed 
subsonic- supersonic flow around this double comer: the outer 

subsonic flow at the sonic boundary will supply the missing 
symmetrical boundary conditions (as for the example in fig. 5(a)). 

But it does not appear impossible that such a mixed suporscnic- 
Bub sonic flow which has been made symmetrical by the outer subsonic 
flow may be unstable. Since the subsonic flow la produced by the 
form of the obstacle in the' flow it will "bo symmetrical only for 
a symmetrical obstacle. On the other hand the outer flow will not 
be symmetrical for an unsymmetrical obstacle, as for instance a 
symmetrical body with a local supersonic region and a compression 
shock. Therefore^ it seams very doubtful whether the outer flow 
whose form is decisively determined by the form of the supersonic 
region, will in turn be able to reshape decisively the supersonic 
flow by forcing the symmetry conditions upon it. 6 

A confirmation of the noneymnetry of the symmetrical type of 
solution would explain the occurrence of unsyrane t ri c al solutions 
with compression shocks which has been observed in tests. 

On this, occasion I should like to point out that a symmetrical 
supersonic flow also oould be suggested for the single wall represented 

6 A solution of this stability problem shall be obtained by 
assuming a small unsymmetrical deformation of the sonic boundary and 
by then observing whether this deformation Increases or diminishes. 
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in figure 8(a) which is physically not real but satisfies the 
conditions for potential flows with regard to continuity and 
irrotatlonallty as well as Bernoulli * b equation. This flow is 
represented in figure 8(h). Such a flow, as is well known, is 
physically real for the reason that a disturbance can not he 
transmitted upstream In supersonic flow. For solutions which 
were obtained analytically it must, therefore, always he especially 
verified whether this additional condition has been satisfied for 
the supersonic flow. 


71. STOMHY 


1. A method of calculating the approximate velocity field 
for compressible flows with local regions of supersonic velocity 
has been presented. Starting from the flow at a large distance 
from the profile determined' according to Pranitl* s rule, this 
outer flow was continued to the sonic velocity boundary by means 
of a numerical method; the method cf characteristics of Prandtl- 
Busemami was applied for continuation beyond that boundary. These 
calculations result finally in the contour of the profile in the 
region of supersonic velocity. 

2. It has been demonstrated in an example that mixed subsonic- 
super sonic flows about two-dimensional bodies can be calculated 
where net only the transition from subsonic to supersonic but also 
the transition from s’personlc to subsonic takes place continuously, 
that is, without pressure Jump. 

3. Ho physically real local area of supersonic velocity could 

be determined for the example considered here when the sound 
velocity was far exceeded because then expansion or compression 
shocks occurred. However, there is a prospect of calculating mixed 
flovB for such cases also: before starting the calculation one 

would have further to extend the outer subsonic flow by means of 
source and sink distributions; in this way the greater need, for 
space of the supersonic velocity field would bo satisfied. 

4. The symmetry of the local area of supersonic velocity for 
mixed supor sonic- subsonic flowB may be enforced by the outer 
subsonic flow. Since, however, the outside subsonic flow in 
turn must be produced by the profile in the flow or by the local 
area of supersonic velocity, the assumption seems Justified that 
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often only an unstable symmetrical flow can be produced by the 
out aide subsonlo flow; and this unstable symmetrical flov vill 
turn even at email disturbances Into the uneymmetrioal case with 
compression shocks. 


Translated by Mary L. Mahl er 
National Advisory Committee 
for Aeronautics 
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Figure 1. Schematic representation of the flow regions for mixed subsonic 
-supersonic flow. 
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Fig. 2 
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Figure 2. Partition of the flow field for the difference method. 
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Sonic Velocity-Boundary 



H) 

Figure 5a. Grid chacteristics in the local supersonic region for a two-dimensional body 
B (relation of thickness d/1 = 0.0715 at the Mach number M = 0.86. 
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Figure 5b. Calculated pressure distribution for two two-dimensional bodies in incompressible 
flow and in compressible flow with local area of supersonic velocity. 
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Figs. 6,7 
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Figure 6. Mixed subsonic-supersonic flow about a body with plane surface' in 
the supersonic region. 



Expansion wave 

Compression wave 


Figure 7. Schematic representation of the two different kinds of wave 
radiation at the sonic velocity boundary. 



Figure 8b. Symmetrical potential flow with supersonic free-stream 

velocity about a physically not realizable double corner. 
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